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Current-driven destabilization of both collinear configurations in asymmetric 

spin-valves 
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Spin transfer torque in spin valves usually destabilizes one of the collinear configurations (either 
parallel or antiparallel) and stabilizes the second one. Apart from this, balance of the spin-transfer 
and damping torques can lead to steady precessional modes. In this letter we show that in some 
asymmetric nanopillars spin current can destabilize both parallel and antiparallel configurations. As 
a result, stationary precessional modes can occur at zero magnetic field. The corresponding phase 
diagram as well as frequencies of the precessional modes have been calculated in the framework of 
macrospin model. The relevant spin transfer torque has been calculated in terms of the macroscopic 
model based on spin diffusion equations. 

PACS numbers: 75.60.Ch,75.70.Cn,75.70.Pa 



Introduction: Spin transfer between conduetion elec- 
tron system and localized magnetic moments gives rise to 
new phenomena, like for instance current-induced switch- 
ing between different magnetic states [l|. Several theo- 
retical models have been developed to describe physical 
mechanisms of the spin transfer and magnetic switching 
phenomena Q, |M |j, [M lyl- Apart from this, magnetic 
switching has been observed in a number of experiments 

At certain conditions electric current can cause tran- 
sition to steady precessional modes, where the energy is 
pumped from conduction electrons to localized magnetic 
moments [H IH El El El 111 . This phenomenon is 
of high importance due to possible applications in mi- 
crowave generation. In typical Co/Cu/Co spin valves, 
the steady precessions exist for external magnetic fields 
larger than the anisotropy field , an d for currents exceed- 
ing certain critical values El Ell- For lower values of 
external field, the current drives switching to antiparallel 
(AP) or parallel (P) states, depending on the initial state 
of the system. 

In this Letter we show that spin transfer torque in 
asymmetric systems (with the two magnetic films having 
different bulk and interface spin asymmetry factors, e.g., 
Co/Cu/Py nanopillars) vanishes at a certain noncoUincar 
configuration due to an inverse spin accumulation in the 
nonmagnetic spacer layer. As a result, spin current can 
destabilize both P and AP magnetic configurations for 
one orientation of the bias current and stabilize both con- 
figurations for the opposite current. The former case is of 
particular interest, as the spin current can excite preces- 
sional modes in the absence of external magnetic field, 
which is of particular importance from the application 
point of view. Such a system is thoroughly studied in 
this Letter. 

Model and description: We consider a tri-layer sys- 
tem which is attached to two external nonmagnetic leads. 
The tri-layer consists of a magnetically fixed (thick) layer 
of thickness di , nonmagnetic spacer layer of thickness d2 , 



and a thin magnetic (free or sensing) layer of thickness 
da . Magnetic dynamics of the sensing layer is described 
by the generalized Landau-Lifshitz-Gilbert equation 
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where s is the unit vector along the spin moment of 
the sensing layer, 7g is the gyromagnetic ratio, /ig is 
the magnetic vacuum permeability, H^s is an effective 
magnetic field acting on the sensing layer, a is the 
damping parameter, and Mg stands for the saturation 
magnetization. Generally, the effective field includes 
an external magnetic field -ffext, the uniaxial magnetic 
anisotropy field iJa, and the demagnetization field Ha] 
iJeff = -^oxt ez - iJa(s • 6^) 6^ -|- iJd(s ' e^) e^, where 
Ex and Bz are the unit vectors along the axes x (normal 
to the layers) and z (in-plane), respectively. The last 
term in Eq.(^) stands for the torque due to spin trans- 
fer, T = T0 -\- T^, with Te = al s X {s X S) = Tgeg, and 
Tip = bl s X S = TipB^. Here, S is the unit vector along 
the spin moment of the fixed magnetic layer {S = e^), 
while eg and e^ are the unit vectors of a coordinate sys- 
tem associated with the polar 9 and azimuthal ip angles 
describing orientation of the vector s. The current / is 
defined as positive when it flows from the sensing layer 
towards the thick magnetic one (opposite to the axis x). 
The current-induced torque r and the corresponding 
parameters a and b have been calculated in the diffusive 
transport regime, assuming that the spin current com- 
ponent perpendicular to the spin moment of the sensing 
layer is entirely absorbed in the interfacial region 2, 3, 6^. 
For the parameter a one then finds jg 
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where g = g//, with the spin accumulation components 
g = idxi gyi dz) taken in the nonmagnetic spacer at the 
very interface with the sensing layer (x -^ Xq), and 
G|| is the mixing conductance |3|. The parameter b is 
given by a similar formula, but with Re{G|j^} replaced 
by -ImjGti} and Im{G^i} replaced by Re{Gti}. 

A qualitative picture of the non-linear dynamics can 
be obtained from the local phase portraits of the corre- 
sponding linearized system in the vicinity of fixed points. 
Accordingly, we linearized Eqs.QJ and calculated eigen- 
values of the corresponding Jacobian. For a two dimen- 
sional problem (as in our case) the eigenvalues depend 
only on the trace and determinant of the Jacobian. Apart 
from this, the eigenvalues are current dependent and de- 
termine stability of the fixed points. More specifically, 
the point is stable when real parts of all the eigenvalues 
are negative, and becomes unstable when at least one of 
them is positive. The system considered, Eqs.(QJ, can 
have several fixed points given by r = figMi^d^s x i?off- 
However, only two of them, P and AP, are trivial with 
position independent of current. In addition, the current 
driven global phase portrait of the system considered be- 
low contains two saddle points with separatriccs dividing 
the phase space to basin of attraction for P and AP fixed 
points and two focuses (located near energy maximum). 

Critical currents: Consider first the P state and as- 
sume positive determinant of the corresponding Jaco- 
bian. Since the eigenvalues of the linearized problem can 
be complex numbers, the non-zero imaginary parts give 
rise to periodic components of the fundamental solutions. 
The condition of vanishing trace determines the critical 
current /™*^ which destabilizes the P state and switches 
system to a precessional (PRC) state (this can be assign 
as Hopf bifurcation with a limit cycle emerging) , 
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where the parameters a and h have to be calculated in the 
limit of ^ 0. With a further increase in current, the 
complex conjugated imaginary parts of the eigenvalues 
vanish at a certain point, at which two real eigenvalue 
branches arise. When current increases further, one of 
the eigenvalues goes to zero, leading to a saddle-node bi- 
furcation. In the case studied below, two new saddles ap- 
pear inside the limit cycle making it unstable [l3|. Thus, 
the critical currents driving system to a static state (SS) 
can be found as, 
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The signs of critical currents depend on the parameters 
a and b taken at — > 0. In the case considered below, 
/ > 0, these parameters obey the conditions a — ba > 
and b < (the following discussion is limited to the 
case, where these conditions arc fulfilled). Consequently, 



the PRC regime holds for I > I^^'^, and the P state is 
stable for / < /™^. Since Jf| < Jf^, the SS states can 
occur for / > I^f_, and we can skip /f^ in the following 
considerations. 

Asymmetric nanopillars: Let us consider now an asym- 
metric nanopillar, where the two ferromagnets and/or 
two external leads arc of different materials, like for in- 
stance Cu/Co/Cu/Py/Cu spin valve. It has been pre- 
dicted recently that the spin-transfer torque acting on 
the sensing layer in such a spin valve can vanish for a 
noncoUincar magnetic configuration [g. In Fig.l(b,c) we 
plot the angular variation of the torques r^ and Te , which 
clearly shows that Te vanishes at 9c ^ 0, tt. Such a behav- 
ior of Te can be observed if an inverse spin accumulation 
builds up in the spacer layer, as shown schematically in 
Fig. 1(a). This can be achieved when the fixed magnetic 
layer has bulk spin asymmetry weaker than that of the 
sensing layer, and its thickness is smaller than the spin 
diffusion length. The general conditions for vanishing 
spin torque at a noncoUinear configuration can be de- 
rived assuming dTe/dO > in both P and AP states, and 
for in-plane configurations {^ -^ '^ 1"^^ it has the form 



±ffz 



Re{GTi} dd 



T 



dgy_ 
do 



<0. 



(5) 



For a general configuration one has to replace gx i— > 
{cjx sin if — (jy cos ip) and ijy i— > (g^ cos ip + gy sin ip). 

The parameter 9c /tt is plotted in Fig.l as a function 
of di and d^ for the right lead made of Cu (d) and Au 
(e), respectively. As one can note, 9c becomes particu- 
larly large for thin Co layers and thick Py ones. How- 
ever, for a too large ratio ds/di the roles of the Co layer 
as a fixed one and Py as a sensing layer could be in- 
terchanged. Therefore, in the following we will consider 
the Cu/Co(10)/Cu(10)/Py(4)/Cu nanopillar (the num- 
bers are the layer thicknesses in nanometers) , with Py as 
a sensing layer and Co as the magnetically fixed one. For 
the Py layer we assume Ms = 10053.1 Oe, i?a = 2.51 Oe, 
Hd = 0.65Ms and a = 0.003. The dependence of 9c on 
the spacer thickness is rather weak, so thicker spacers can 
be used to eliminate possible intcrlayer exchange interac- 
tion. From comparison of Figs 1(d) and 1(e) follows that 
reducing the spin-flip length in the lead adjacent to the 
sensing layer increases 9c and also enhances the torque 
acting on the free layer (reducing the relevant critical cur- 
rent). This is due to an increased spin accumulation at 
the Cu/Py interface. The t^ component is usually much 
smaller than Te, and therefore plays a negligible role in 
the initial stage of the switching from P state. 

From the above follows that in asymmetric spin valves 
both P and AP states can be unstable when current den- 
sity exceeds a certain critical value. Accordingly, preces- 
sional states and/or bistability of the static states can be 
expected for certain values of spin current [ij [13, llal ■ 
If energy of a static state is close to the magnetic en- 
ergy maximum, and a critical current for its stabiliza- 
tion is larger than the critical current needed to destabi- 
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FIG. 1: (color online) (a) Schematic structure of the nanopil- 
lar and spin accumulation in thin spacer layer for the par- 
allel configuration of the magnetic layers, and for elec- 
trons flowing along the axis x. Angular dependence of 
the torques t^ (b) and re (c), acting on the Py layer in 
the Cu/Co(10)/Cu(10)/Py(4)/Cu spin valve (independent of 
if), (d) 9c /n as a function of di and da in the nanopil- 
lar Cu/Co(di)/Cu(10)/Py(d3)/Cu. The white dashed line 
corresponds to a critical thicknesses for ZFOs (zero-field os- 
cillations), (e) 6»c/7r in the Cu/Co(di)/Cu(10)/Py(d3)/Au 
structure. The other parameters are: (Co) bulk resistiv- 
ity p* — 5.1/if2cm, spin asymmetry factor f} = 0.51, and 
spin-flip length Isi — 60 nm; (Py) p* = 16/^r2cm, /? = 0.77, 
IbI = 5.5 nm; (Cu) p* = 0.5 fj.Qcm, Ibi = 1000 nm; (Au) p* — 
2 /ificm and lat = 60 nm. In turn, for the Co/Cu interfaces we 
assume the interfacial resistance per unit square R* — 0.52 x 
10~^'''Slm^, interface spin asymmetry factor 7 = 0.76, and 
the mixing conductances Re{Gtx} = 0.542 x 10^^ il~^m~^, 
and lm{Gn} = 0.016 x lO^^fi^^m"^; for the Py/Cu inter- 
faces we assume R* = 0.5 x 10"^^ fim^, 7 = 0.7, RelGn} = 
0.39 x 10^^ n-^m-^ and Im{Gu} = 0.012 x lO" H" V'^ 



lizc the state corresponding to the energy minimum, the 
PRC state is allowed |Tl|. Thus, the sufficient condition 
for the presence of oscillations reads /f^ > /(f^^. The 
above condition, together with the necessary condition 
of vanishing spin torque in a noncoUinear configuration, 
allow us to determine the region where the zero-field os- 
cillations (ZFOs) exist. In Fig.l(d-e), the white dashed 
lines (corresponding to /f^ = I^^'") separate the regions 
where the ZFOs occur from those where they cannot ex- 
ist. We note that critical currents for switching from the 
AP state are still given by Eqs lO and l@}, but with r.h.s. 
multiplied by —1 and Hcxt replaced by — i/cxt. 

Figure 2 shows dynamical phase diagram of the 
Cu/Co/Cu/Py/Cu nanopillar. The reduced magnetore- 
sistance, r = {1 - cos^(6'/2))/(l -I- cos^(6'/2)), is shown 
there as a function of the external field and reduced cur- 
rent density I/Iq, with Iq — lO^Acm"^. The dashed 
lines correspond to the critical currents given by the an- 



alytical formulas ^ and Q , and fit very well to the nu- 
merical results. For weak external fields we find a stable 
static state with a high magnetoresistance value (HSS), 
which is close to maximal magnetic energy. The mi- 




FIG. 2: (color online) Dynamical phase diagram of the magne- 
toresistance for the Cu/Co(10)/Cu(10)/Py(4)/Cu spin valve, 
scanned by increasing (left) and decreasing (right) current 
with the sweeping rate 59 A/cm^s in a constant applied field. 
The dashed lines correspond to critical currents given by Eqs 
@ and 1^. The other parameters as in Fig.l. 




FIG. 3: Fundamental frequency of the magnetoresistance 
oscillations in Cu/Co(10)/Cu(10)/Py(4)/Cu, calculated as 
a function of the reduced current swept with the rate 16 
A/cm^s in a constant external field as indicated. The in- 
set shows precession orbits in zero external magnetic field for 
J/Jo — 0.45, 0.5, 0.6, 0.7. The other parameters as in Fig.l. 



crowave oscillations exist in the area between the two 
dashed white lines in Fig. 2. (left part). As one can see, 
these oscillations exist in the absence of magnetic field 
for a certain current region. The fundamental frequency 
ujo (first harmonics) of the magnetoresistance oscillations 
is shown in Fig. 3 as a function of the current density and 
for different magnetic fields. Along the segments, where s 
moves almost in the layer plane, the spin processes mainly 
around _ffa and -ffcxt, with the angular velocity propor- 
tional to \jg\{Ha + Hcxt), see the inset in Fig. 3. Along the 
remaining part (s moves almost perpendicularly to the 
layer plane), the angular velocity is larger and propor- 
tional to |7g|(i:fa + ^cxt + -f^d)- With increasing /, the 
average orbital speed decreases while the arc length of 



the orbit increases. Consequently, luq decreases with in- 
creasing /. Such precessions of decreasing frequency with 
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FIG. 4: Current driven hysteretic behavior of the magnetore- 
sistance and spin components of the sensing layer in zero field, 
scanned with increasing (solid line) and decreasing (dashed 
line) current with the sweeping rate 140 A/cm^s. The pa- 
rameters as in Fig. 2. 

increasing current in symmetrical spin valves are called 
'clamshell' [13 or 'in-plane' |16J modes. But the main dif- 
ference between symmetrical pillars and those studied in 
this Letter concerns the way the spin transfer acts on 
s. In our case, increasing current leads to an increase 
in X, and the orbits become narrower and elongated in 
the z direction due to Tg. At the critical current den- 
sity given by Eq.(@J), when Sz — cos 0c, Te vanishes and 
orbital velocity reaches minimum. The damping torque 
is small, and the orbit bifurcates into SS states (fixed 
points) due to a non-zero r^. A further current increase 



at low fields drives the system via a transient regime to 
HSS state with positive or negative s^ component (see 
Fig. 4). The presence of HSS states gives rise to a current- 
driven hysteresis (Figs. 2 and 4). The hysteresis for lower 
fields can be formally written as a sequence of transi- 
tions P^PRC^SS^T+^HSS^T_^SS for increasing 
current and SS^T+^HSS^T_^OOP->P/AP for the 
decreasing current, with final bistability reached via the 
'out-of-plane' precessions (OOP), observed also in sym- 
metrical spin- valves ^^. The switching from SS state 
to a HSS state is mediated by the transient regime T_|_ 
(where spin transfer pumps energy to the system). The 
HSS states can be observed in a quite broad region, up 
to I /Iq ~ 22.5; and further current increase switches sys- 
tem via the next transient T_ regime (where spin transfer 
dissipates energy) to the SS state. For higher fields the 
HSS states become unstable and the hysteresis vanishes. 
Finally we note that for r^p = the static states will not 
occur and the area of PRC states extends for / > /™'^ . 
In conclusion, we have shown that spin polarized cur- 
rent in certain asymmetrical structures can destabilize 
both P and AP configurations of a nanopillar spin- valve, 
and drive microwave oscillations in the absence of mag- 
netic field. This takes place only for one orientation of 
the bias current. We have also presented the correspond- 
ing dynamical phase diagram. 
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